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Abstract 

In [3] A.D. Aleksandrov raised a general question of finding variational statements 
and proofs of existence of polytopes with given geometric data. The first goal of this 
paper is to give a variational solution to the problem of existence and uniqueness 
of a closed convex hypersurface in Euclidean space with prescribed integral Gauss 
curvature. Our solution includes the case of a convex polytope. This problem was also 
first considered by Aleksandrov and below it is referred to as Aleksandrov's problem. 
The second goal of this paper is to show that in variational form the Aleksandrov 
problem is closely connected with the theory of optimal mass transport on a sphere 
with cost function and constraints arising naturally from geometric considerations. 
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1 Introduction 



In his book on convex polytopes A.D. Aleksandrov raised a general question of finding 
variational statements and proofs of existence of polytopes with given geometric data [3], 
section 7.2.4. As examples of geometric problems in convexity theory for which variational 
solutions are possible Aleksandrov referred to the Minkowski problem for which such a proof 
was given by Minkowski himself [25] (see, also, Schneider [34], section 7.1 and Klain [18]) and 
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the Weyl problem for which a variational approach was sketched by Blaschke and Herglotz 
[5] and later carried out by Volkov [37]. 

The first goal of this paper is to provide a variational solution to the problem of 
existence and uniqueness of a closed convex hypersurface in Euclidean space with prescribed 
integral Gauss curvature. Our solution includes the case of a convex polytope. This prob- 
lem was first considered also by Aleksandrov and we refer to it as Aleksandrov's problem. 
Aleksandrov studied it with a nonvariational approach in [1] (see also [3], ch. 9) using his so 
called "mapping lemma" which is a variant of the domain invariance theorem and requires 
uniqueness of the solution in order to establish existence. By contrast, our solution of Alek- 
sandrov's problem, in addition to being variational, does not require uniqueness for the proof 
of existence. 

The second goal of this paper is to show that in variational form Aleksandrov's prob- 
lem can be considered as a problem of optimal mass transport on a sphere with appropri- 
ate cost function and constraints. As a result of this connection with optimal transport 
"economics-like" interpretations of some classical geometric concepts such as the Gauss map 
are obtained. Furthermore, this connection permits, in principle, a numerical determination 
of a convex polytope by the integral Gauss curvature by methods of linear programming. 
This should be useful in some applied problems involving numerical determination of convex 
surfaces by Gauss curvature [20, 19, 33, 28]. 

In order to state the Aleksandrov problem and outline our main results we recall 
the notion of integral Gauss curvature. In Euclidean space R n+1 , n > 1, fix a Cartesian 
coordinate system with the origin at some point O. By a closed convex hypersurface we 
understand here the boundary of a compact convex set in R n+1 containing O in its interior. 
The set of such closed convex hypersurfaces is denoted by T n . Let S n be a sphere of unit 
radius centered at O and F e T n . Since F is star-shaped with respect to O, we can 
parametrize it as r(x) = p(x)x, x e S n , where p(x) is the radial function giving the distance 
from O to the point of intersection of F with the ray of direction x originating at O. (Here 
and elsewhere in the paper a point on S ra is treated also as a unit vector in R n+1 originating 
at O.) The Gauss map Np : F — > S n maps a point X e F into Np(X) C S n , which is the 
set of outward unit normals to all supporting hyperplanes to F at X. Define the generalized 
Gauss map as a F = N F o r : S n — > S n . At x G S™ such that at r(x) G F there is more than 
one supporting hyperplane the map dp is multivalued. If u is a subset of S n then we set 
ap{uj) = \J xeu> a F{ x )- For any Borel subset wcS" the set a F (uj) is Lebesgue measurable on 
S n [34], section 2.2. The function a(a F ), where a is the standard n— dimensional Lebesgue 
measure on S n , is a finite, nonnegative and countably additive measure on Borel subsets of 
S n . It is called the integral Gauss curvature of F (transferred to S n ) [1]. 

The Aleksandrov problem posed by him in [1] is to find conditions on a given measure 
H on S n under which there exists an F G T n such that 

a{a F {uo)) = fx(u) for any Borel set to C S n . (1) 
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In the same paper Aleksandrov proved the following 

Theorem 1. In order for a given function /i on Borel subsets ofS n to be the integral Gauss 
curvature of a convex hypersurface F G T n it is necessary and sufficient that 

(i) ji is nonnegative and countably additive on Borel subsets of S n ; 

(ii) /i(S") = <x(S")/ 

(iii) the inequality 

li(9 n \u)>a(u*), (2) 

holds for any spherically convex subset uj C S n , uj ^ S"; here, uj* is the dual to uj, that is, 
u* = { y e S" I (x,y) < G uj}. 

The necessity of (i)-(iii) is easy to verify. The proof of sufficiency given by Aleksandrov 
consists of two steps. In step one the theorem is established for convex polytopes when the 
function p is an atomic measure concentrated at a finite number of points x±, ...,xk G S n . 
Here, to prove existence Aleksandrov uses his mapping lemma proving first uniqueness among 
polytopes up to a homothety with respect to O. In step two the initially given function p is 
approximated weakly by a sequence of atomic measures for which, by step one, there exists a 
sequence of polytopes in T n solving the problem in this case. It is shown that this sequence 
contains a subsequence converging to a convex hypersurface in the Hausdorff metric, the 
limiting hypersurface is in T n and its integral Gauss curvature is p. Uniqueness in J rn up to 
a homothety with respect to O was proved by Aleksandrov in [2]. Later, Pogorelov [30], [31], 
[32] studied the Aleksandrov problem and its generalizations by a different approach not 
relying on the mapping lemma. The variational solution presented below is different from 
that by Pogorelov and gives significantly more information on the solution and its properties. 

In this paper the hypersurface required in Theorem 1 is obtained as a minimizer for 
the functional 

Q[h,p]= [ logh(N)da(N) - [ \ogp{x)dp{x) (3) 

on the set 

A={(h,p)e C(S n ) x C(S n ) | h > 0, p > 0, (4) 
log h(N) - log p(x) > log (x, AT), Or, iV) G S n x S™, (x,N) > 0}. 

It is shown that a minimizing pair (h, p) exists, unique (up to a multiplicative constant) 
and h and p are, respectively, the support and radial functions of a hypersurface F G T n 
satisfying (1). In agreement with the question of Aleksandrov asking for a variational solution 
in the class of convex polytopes we first solve the minimization problem for polytopes and 
then prove existence of minimizers in the general case. As can be seen from considerations in 
section 2, the choice of the functional Q and of the set A is a completely natural consequence 
of the geometric duality relating the radial and support functions of a convex hypersurface 
in T n via a Legendre-like transform. 
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It turns out that the generalized Gauss map ap of the minimizer of Q maximizes the 
functional 

/ c(0- 1 (N),N)da(N), 9eQ, (5) 
where the "cost" function is given for (x, N) G S n x S n by 

c(x N) = l lo § W> when (x, iV) > 
^ ' ' [—00 otherwise ^ ' 

and 6 = {9} is the set of measurable maps of S n onto itself, possibly multivalued, satisfying 
the following conditions: 

the image 9(u) of any Borel set u C S n is Lebesgue measurable, (7) 
a{N G S ra I 9~ 1 (N) contains more than one point} = 0, (8) 

and for fi as in Theorem 1 and any / G C(S n ) 

/ f(9-\N))da(N) = [ f(x)d^x). (9) 
Js n Js n 

In addition, the maximizer ap is unique (up to sets of measure zero with respect to //). 

From the geometric point of view, the function \og(x, Np(x)) gives a scale invariant 
quantitative measure of "asphericity" of a hypersurface F with respect to O. For example, 
for a sphere centered at O it is identically zero, while for a sufficiently elongated ellipsoid 
of revolution centered at O it has large negative values at points where the radial direction 
is nearly orthogonal to the normal. Roughly speaking, the above result says that the most 
efficient way (with respect to the cost c(x, N)) to transfer to a an abstractly given measure \i 
on S" is to move it by the least possible distance and this is accomplished by the generalized 
Gauss map ap of the convex hypersurface F solving the problem (3), (4). 

The above result establishing that ap maximizes the functional in (5) can also be 
viewed as a counterpart of a result by Brenier [6], [7] (and for more general cost functions by 
Gangbo-McCann [10], [11] and Caffarelli [8]) showing that among measure preserving maps 
between two sets U and V in Euclidean space the map optimal with respect to the cost 
function \x — y\ 2 , x G U, y G V, is unique and is the subgradient of a convex function. 

In the framework of mass transport theory on S n the problem (3), (4) is the dual of 
the following primal maximization problem of finding 7 such that 

C[y] = sup C[y], where £[7] := / / c(x, N)d>y(x, JV), (10) 

r(/x,<r) JS" Js n 

and r(/x, a) is a set of joint Borel measures on S"xS n with marginals \x and a as in Theorem I. 
In this case, the optimal measure is generated by the map ap, that is, ^[U, V] = a[ap(U)nV] 
for any Borel subsets U and V on S n . In addition, the duality relation 
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holds. Thus, similar to the classical theory of optimal mass transport developed by L.V. 
Kantorovich [16], [17] and many other authors, we also have here the primal and dual 
problems and the usual duality between them. In this respect, the problem (10) can be 
viewed as a direct generalization of the classical problem. On the other hand, the problem 
(10) is essentially different from the classical problem. In our case the transport takes place 
on a sphere and the cost function (6) does not have a direct analogue among cost functions 
considered previously in Euclidean space [10], [11], [8] or in Riemannian space [9]. Cost 
functions of the type (6) seem to be more common (and natural) in geometric problems than 
in economics or mechanics. 

It is worthwhile noting that the variational approach used in this paper is quite general 
and is expected to be useful for studying other nonlinear geometric problems, in particular, 
the generalized Minkowski problem considered in [21], [22], [23]. In a simpler situation 
this approach was used in [29], and in [13], [14], [12] it was applied to study problems in 
geometrical optics. 

Concerning regularity of the solution to the Aleksandrov problem, we note that in 
analytic form, when the measure \i has a density m{x) > 0, x G S n , the problem reduces to 
finding the radial function p satisfying the following equation of Monge-Ampere type [27] 

(P 2 + I Vpir^P 1 - ^'^" + 2p f J + ^ =mon S«. (11) 

Here, the gradient V is computed in the standard metric e = ^2™j =1 eijdu l dui on S ra , u 1 , ...,u n 
are some local coordinates on S n , Vijp, i,j = l,...,n, are the second covariant derivatives 
in metric e, while pi = §^,i = l,...,n. Existence of smooth solutions to the Aleksandrov 
problem under appropriate smoothness assumptions was established for n = 2 by Pogorelov 
[31], Chapter VIII, and for arbitrary n by the author [27]. For related results concerning C° 
apriori bounds and stability of solutions see Kagan [15] and Treibergs [35]. An analogue of 
Aleksandrov's problem for nonparametric convex hypersurfaces defined over R n was studied 
by Aleksandrov in [2], Pogorelov [30, 31, 32], Bakelman [4] and, more recently, by McCann 
[24] and the author [29]. 

The rest of the paper is organized as follows. In section 2 we study a Legendre- 
like transform between radial and support functions of hypersurfaces in jF n , introduce the 
class of admissible functions and formulate the minimization problem for the functional Q. 
In section 3 we specialize this minimization problem to the class of polytopes and prove 
existence and uniqueness (up to rescaling) of minimizers. In section 4 we establish existence 
and uniqueness of a minimizer to the minimization problem for Q in the class of general 
closed convex hypersurfaces. Finally, in section 5 we describe the connections between the 
variational solution of Aleksandrov's problem and the mass transport theory. 

The author thanks the reviewer for careful reading of the manuscript and many helpful 
comments. 
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2 A variational formulation of Aleksandrov's problem 



The key step leading to a variational statement of the problem is based on the simple but 
crucial fact that two representations are available for any F G T n and a transformation 
generalizing the classical Legendre transform connects these two representations. Let F G 
jF n , p its radial function and r(x) = p(x)x, x G S n , its position vector. An alternative 
representation of F can be given in terms of its support function defined as 

h(N) = sup p(x)(x,N), N G S n . (12) 
xes n 

Since O G intB(F), where B(F) is the convex body bounded by F, functions p and h are 
both positive on S n . It follows from (12) that 

wrsiwj'^ 8 "- (13) 

Thus, the functions p and h provide two descriptions of the same hypersurface F and these 
two descriptions are related by (12) and (13). We refer to (12) as the generalized Legendre- 
Fenchel (LF) transform of p and write h = LF(p). Similarly, the function p defined by (13) 
is considered as LF~ l (h). 

Theorem 2. Let p > and h > be two continuous functions on S n satisfying (12) and 
(13). Then there exists a unique closed convex hypersurface F G T n with radial function p 
and support function h. 



Proof. Define the function 

H(u) = sup p(x)(x,u), u G R n+1 . 

Obviously, H{u) = h(u) for u G S n and, thus, H is a positively homogeneous extension of 
order one of the function h from S n to R n+1 . It is clear that H is also subadditive. By 
a theorem of H. Minkowski it is the support function of a unique compact convex body in 
pri+i. see [34] ? section 1.7. The boundary of that body is a closed convex hypersurface F. 
Since H{u) > for u ^ 0, it is clear that F G T n . 

Let 

R(v) = sup{A >0\ \ve B(F)} for v G R n+1 \ {£>} 

be the radial function of the convex body B(F) (see [34], section 1.7.) Then R is positively 
homogeneous of order —1 and it is defined by its values on S n . Furthermore, R(x)x G F for 
any x G S n . 

We show now that R(x) = p(x) for all x G S n . Fix some N G S n and consider the set 
G N = {x G S n | R(x)(x,N) = h{N)}. 
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Note that G N ^ 0. Since R, h > on S n , (x, N) > for all x G G N . Then by (12) R(x) > 
p(x) for all x G G^. Suppose there exist x G Gn such that p(x) < R(x). Since R(x)x G F, 
the point p(x)x is an interior point of the convex body B(F). Then p(x)(x,u) < h{u) for 
any u G S n and, consequently, 

1 (x, -u) 

"T^r > sup — — . 
p(x) ueS ™ h(u) 

But this contradicts (13). Thus, p(x) = R(x) for all x G Gat. Since U^ves™ = S n , we 
conclude that R(x) = p(x) for all x G S n . QED. 



In contrast to the classical Legendre transform (conjugation) defined by 

/*(«) = sup{( M ,a;) - /(x) I x G R n }, u G R n , 

where / : R™ — > R is a closed convex function (see [34], section 1.6), the LF transform 
relates p and h multiplicatively which is inconvenient for defining the required functional in 
the Aleksandrov problem. It is more convenient to work with the logarithms of functions p, 
h and (x, N), which are all positive for F G T n . Put 

p := logp, h := log/i, K{x,N) := log (x, N) for x, N G S" and (x, N) > 0. 

Then, it follows from (13), (12) that for a F G T n the generalized Gauss map can be defined 

as 

a F (x) = {N G S" I (x, N}>0 and ft(iV) - p(x) = K(x, N)}, x G S n . (14) 

Similarly, 

ap\N) = {x G S" I (x, iV) > and h(N) - p{x) = K(x, N)}, N G S n . (15) 



In this notation the set of pairs of admissible functions (4) is 

A = {(h,p) G C(S n ) x C(S n ) \h>0, p > 0, (16) 
h(N) - p(x) > fc(x, N), (x, N) G S" x S n , (x, N) > 0} 

and the functional (3) is 

Q[h,p}= [ h(N)da(N)- [ p(x)dp(x), (h,p)eA, (17) 

where p is as in Theorem 1. The variational problem we now consider is to determine a pair 
(h, p) G A such that 

Q[M = infg[M- (18) 

A 
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3 A variational solution of Aleksandrov's problem for 
polytopes 



Just as in the classical case, it is convenient to investigate the variational problem above 
first in the class of convex polytopes in T n . Let X K := {xi, ...,xk} be a set of points on 
S n not contained in one closed hemisphere. In this section the set X K will remain fixed. 
Let h, ...,1k denote the set of rays originating at O and going through the points xi, xk- 
We denote by T K the subset of T n consisting of the convex polytopes with vertices possible 
only on the rays Ii,...,Ik- Let P G T K . Denote by T\,...,Tk the points of intersection of 
h, ...,1k with P. It is useful to recall that P is the boundary of the convex hull of the points 
7*1, ...,rK. Thus, for the polytope P its radial function p is completely defined by the vector 
(p(xi), ...,p(x K )). Of course, r» = p{xi)Xi, i = 1, K. 

For the rest of the paper, when speaking of a convex polytope in T K with radial 
function p and if there is no danger of confusion, we put for brevity pi := p(xi), i — 1, K. 

Note also that for a Borel subset wcS" the integral Gauss curvature a(ap(u)) = 
if Xi uj for any i — 1, K. 

We will need a version of Theorem 2 for polytopes. 
Proposition 3. Let X K be as above, pi > 0, ...,Pk > 0, and h G C(S n ), h > 0. Suppose 

h(N) = sup Pi(xi,N), NeS n , (19) 

l<i<K 

~= SU P T7W' < = l>->*- (20) 
Pi wes™ h(N) 

Then there exists a unique convex polytope P G T K with support function h and radial 
function p such that pi = Pi, i = 1, K . 



Proof. Consider the convex hull P of the points p\X\, ...,pkXk- Because x±,...,xn are 
not contained in any closed hemisphere of S™ and pi > for alH = 1,...,K the polytope 
P G T K . Let pp(x), x G S n , and h P (N), N G S™, be, respectively, the radial and support 
functions of P. Then pp and hp satisfy (12) and (13). It follows from (12) and (19) that 
hp(N) > h{N) VW G S n . On the other hand, for any fixed JVeS" and any x G S n such 
that pp(x)(x, N) = h P (N) we have 

h P (N) = N) < h(N), \ t > 0, J2 X i = !> 

where the sums are taken over the vertices of P whose linear combination gives p P (x)x. 
Thus, hp(N) = h(N) ViV G S n . The equalities pp(xi) = Pi for % — 1, ...,K are proved using 
arguments similar to those at the end of the proof of Theorem 2. QED. 
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For a vector p = {pi, ...,Pk) £ R K , we write p > when > 0, ...,Pk > 0. Also, for 
p > put p = {pi, ...,p K ), where pi := logpj, i — 1, X. Define now the admissible set as 

^ = {(h,p) G C{S n ) x | h > 0, p> 0, (21) 
h{N) -p t > £( Xi , A), x t e X K , A G S n , ( Xi , N) > 0, % = 1, K}. 

Put 

K 

Q K [h,p] = / h{N)da{N) - Vp./i,, G ^, (22) 

where fii, ...,[Ik are positive numbers. When P is a convex polytope in J-"£ with the radial 
function p, we write QK[h,p] setting in (22) p = {pi, p K ). 

Theorem 4. Let Pi> 0, i = 1, K, 

K 

5> = a(S") (23) 
i=i 

and for any (solid) polyhedral convex cone V ^ R n+1 with vertex at O, possibly degenerate, 

J2»i>°((S n nV)*), (24) 

where the sum is taken over all i such that X{ G" S n PI V . Then there exists a polytope P G 
with radial function p and support function h such that 

Q K [h,p\ = inf Q K [h,p\ = inf Q K [h,p\. (25) 
The pair (h,p) is unique up to rescaling {ch, cp) with any c = const > 0. 



Proof. First, we note that it suffices to look for a minimizer of Qk over J-"£. Indeed, if 
{h,p) G Ak is such that 

h(N)-p i >Kl(x i ,N) 

for some i and all JV G S n such that {xi, N) > then we can increase pi slightly to p\ and 
decrease h to h' for some iV so that the pair {h',p') is still in Ak- The functional Q K will 
not increase under such change. Therefore, it suffices to minimize Qk only on pairs in Ak 
such that 

h{N) = sup\pi + iC{ Xi , A% AG S n , (26) 

i 

-Pi = sup [-h{N) + JC{xi, A)], i = 1, K, (27) 
iVeS" 

with the supremum achieved in (26) for each A at some i and in (27) for each i for some 
A. By Proposition 3 each such pair {h,p) defines a unique convex polytope in T K with the 
support function h and radial function p such that pi = Pi, i = 1, K. 
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It follows from (23) that if a pair (h,p) G Ak and c is a positive constant, then 



Q K [ch,cp] = Q K [h,p]. 



Hence, the search for a minimizer in T K can be restricted to convex polytopes in T K such that 
Pi < 1 for alH G {1, K} and pj = 1 at least for one % G {1, K}. To avoid introduction of 
additional notation, we will continue to denote by T K the set of convex polytopes satisfying 
this restriction, while T K will denote the closure of the set of such polytopes with respect 
to the metric 



dist(P l ,P 2 ) 



K 



£(p!-p*) 2 



i=i 



1/2 



The set of vectors {(pi, Pk)} corresponding to polytopes in T K forms a closed bounded set 
in R^ while the corresponding support functions form a compact set in C(S n ) with repect 
to the uniform convergence on S n . Note that the functional Qk is continuous on T K - 

Remark 5. Since the set Xk is not contained in a closed hemisphere, it is easy to construct 
a polytope in T K on which 

\Q K [h,p}\ <oo. (28) 
For example, the polytope inscribed in S n with vertices xi, ...,xk satisfies (28). 



Next, we show that the inf^n Qx[h, p] is attained on a closed convex polytope for 
which pi > Vi = 1,...,K. Since X K is not contained in a closed hemisphere, on this 
polytope h(N) > ViV G S n and then Qx[h, p] > — oo. Consequently, the Remark 5 and 
the continuity of Qk on T K imply the existence statement of the theorem. 

Suppose the inf^n Q.K[h,p] is attained on a convex polytope P = (h,p) and O G P. 
Again, because Xk is not contained in a closed hemisphere, at least one vertex of P must be 
at O. The radial functions of the polytopes in T K are uniformly bounded and there exists 
a sequence of polytopes P e = (h e ,p t ) G T K converging to P as e — > 0. Note that for any 
e < 1 the inequality p\ < e cannot be true for all % — 1, ...,K as p\ — 1 for some i. The 
same holds also for p 1 , ...,p~K- Let {x 1 , ...,x s } be the set of directions for which p\ — > and 
{x s+ i, ...,xk} the set of directions for which p\ y4 0. We have 

Q K [h\ P e ]= [ V(N)da(N)-J2pi»i = J2\ f [p\ + £{xi,N)]da{N)-p\p\ (29) 

s K K „ 

= J]/ 3 i^(«P«(^))-^]+ A e [^(«^(^))-^] + ^ / t(xi,N)da(N). (30) 

i=l i=s+l i=l ^f 1 ') 

We consider now separately each of the terms in (30), beginning with the first term on the 
left. 

Any polytope P e is the boundary of the convex hull of vertices of P e on the rays 
l s +i, ...,1k and of the polytope P{ which is the boundary of the convex hull of O and vertices 
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of P e on the rays l±, l s . As e — > the polytopes Pf contract to O while P e converges to 
P. Denote by V the polyhedral angle formed by the faces of P adjacent to O and observe 
that V is the boundary of the convex hull of the rays Z s+ i, ...,1k- In addition, 

s 

J2cr(a pc (x t )) - = a((S n fl V)*). 

i=i 

It follows from (24) that for e sufficiently small we must have 

s s 

^to> ^2(j(ape(xi)). 
i=i i=i 

On the other hand, for any i = 1, s the convex polytopes with vertices {rj = p e jXj, j = 
i, s + 1, fT} also converge to P. Consequently, 

s s 

lim n Yl Pi [v{up*{xi)) ~ to] = !im( min pf) V [a(a P e( Xi )) - p t ] = +oo. (31) 

1=1 i=l 



Since p- > 5 > for some 5 and all e and alii = s + 1, K, for the second sum in 
(30) we have 

K 

PiW( a A x i)) ~ to 

i=s+l 



< oo for all e. 



(32) 



Finally, we consider the last term in (30) and estimate its absolute value from above 
by a bound independent of e. Observe that for each P e and any Xi G X K the inequality 
(xi,N) > holds for all N e ape( Xi ), since the origin (9 is strictly inside the convex body 
bounded by P e . Using spherical coordinates on S n with the origin at Xi, we obtain 



/ 

J a 



apt (xi) 



\fc(xi,N)\da(N) 



{xi,N)>0 



where Xa pt (x l ) is the characteristic function of the set ape( Xi ) and v n -i is the (n — 1)- 
dimensional volume of S™" 1 . Then 



E 



t{xi, N)da(N) 



1 Ja P e(xi) 



< 



Kv n _\ when n > 2 
Kn/2 when n = 1. 



(33) 



It follows from (29)-(33) and remark 5 that inf^-n Qx[h, p] is attained on a polytope 
in JF]^. This completes the proof of the existence statement in this theorem. We now prove 
uniqueness of the minimizer of Qk (up to a homothety with respect to O). 
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Let P 1 = (h l , p 1 ) and P 2 = (h 2 ,p 2 ) be two minimizers of Qk in T K . After rescaling, 
if needed, we may assume that p 2 < pj, Vi = 1, K with the equality holding at least for 
one i and strict inequality for some other i. Then 

Q K [h\ p 1 ] < Q K [h\ p 2 } 

and we arrived at a contradiction. QED. 

Theorem 6. Let P be a minimizer of Qk in Theorem 4- Then 

a(ap(xi)) = pii Vi = 1, ...,K. (34) 

Proof. Suppose that (34) is not true for some m. We show that in this case the value of 
the functional Qk can be reduced. 

Note first that because of (23) and because 

K 

^2a(ap( Xi ))=a(S n ) VP G T K 

i=l 

it may be assumed that cr(ap(x m )) > p m . 

Consider a convex polytope P s which is the boundary of the convex hull of points 
Ti = pfxi, 1 < i < K, where pf = pi when i ^ m and p 5 m = p m — 5, and 5 > is sufficiently 
small so that r, G P s Wi and P s G T K . This is possible since a(ap(x m )) > p m > 0. (Here 
the restriction Pi < 1 Vrcj G X K with the equality achieved at least for one i imposed on 
polytopes in T K in the proof of Theorem 4 is not imposed.) The support function of P s is 
denoted by h 5 . 

Clearly, 

K K 



{Jap(xi) = {Japs(xi) = S n 



i=l i=l 



and a P s(x m ) C ap(x m ), while ap(xi) C a P s(xi) for i ^ m. Put = a P s(xi) fl ap(xj). 
Then 



aps(x m ) = ap(x m ) \ I |J K im J , a P s(xi) = ap(xi) U K im for i 7^ 



m. 



For i ^ m and iV G ap(xj) we have 



h\N) = sup [p* + K(xj, N)] =p d t + K{x u N) = sup [~ Pj + K(xj, N)] = h(N). 

l<j<K l<j<K 

12 



For N e a P s(x m ) we have 

h s (N) - h(N) = sup [p* + £( Xj , N)} - sup [Pj + fc( X j, N)\ 

l<j<K l<j<K 

= pL + £(x m , N) - sup [pj + K(xj, N)] < p s m -p m = -— + o(5). 

l<j<K Pra 

In addition, for suitable i (depending on N) 

h s (N) - h(N) =pf + iC(xi, N) - sup [pj + fc( Xj , N)\ < ViV e re™, i ^ m. 

Then 

/ \h 5 (N) - h(N)]da(N) = [ [h s (N) - h(N)]da(N) (35) 
+ [ $ S ( N ) - HN)}da(N) = f [h 5 (N) - ~h(N)]da(N) 
+ ih S (N) - h(N)]da(N) < (~ + o(5)) a(a pS (x m )). 

i+m J«.im V Pm J 

It is shown in [15] that 

a(ap(x m )) = a(a P s(x m )) + C5 + o(8), 
where the constant C depends only on the polytope P. Using this in (35) gives 

/ [h 5 (N) - K(N)]da(N) < —L a (a P (x m )) + o(5). 

Then for sufficiently small 5 > 

~ 5 

Q K [h 5 , p 5 ] - Q K [h, p] < -— K«p(x m )) - p m ] + o(5) < 0. 

Pm 

However, this is impossible since inf^x. Q K = Qx[h, p\. QED. 



4 The general case 

We use now the above results for convex polytopes to solve the variational problem (18). 
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Theorem 7. Let fj, be a measure on S n satisfying conditions (i)-(iii) in Theorem 1. Then 
there exists a closed convex hypersurface F G T n with support function h and radial function 
p such that 

Q[h,p}=m(Q[h,~p]- (36) 
In addition, for any Borel set wcS" 

a(ap(uj)) = (37) 

that is, the minimizer F is a solution of the Aleksandrov problem. 

Proof. Partition S n so that S n = \Jf =1 V u where V* are Borel subsets of S n such that 
VidVj = when i ^ j and diam(Vi) < 5. Pick a point Xi G V*, and put pi = p{Vi). It is shown 
by Aleksandrov in [1] that for S sufficiently small the set Xk = • i/c} and the numbers 
/ii, ...,Hk (or a subset of X K corresponding to positive /ij) satisfy the conditions of Theorem 
4. Therefore, for each sufficiently large K there exists a polytope P K G J 7 ^ with support 
function h K and radial function p K such that Q K [h K , p K ] < Q K [h K , p K ] V(h K , p K ) G 
By rescaling, if needed, it can be assumed that the diameters of all P K are equal 1. Therefore, 
for 5 — > the corresponding sequence {P K } contains a converging subsequence, which we 
denote again by {P K }- The limit of this subsequence we denote by F and its support and 
radial functions respectively by h and p. 

It is also shown in [1] that the origin O is strictly inside the convex body bounded by 
P and as 5 — > the integral Gauss curvatures of P K converge weakly to the integral Gauss 
curvature of F. This implies (37). The above also implies that 

Q K \h K ,p K ]^Q[h,p}- 



It is known ([34], section 2.2) that for any F G T n 

a{N G S n | ap\N) is multivalued} = 0. 

It follows from (34) that, in particular, 

k „ K 

E / f{af K {N))da{N) = ^ f(x)d» K (x), V/ G C(S n ), (38) 



i=l "<XpK (xi) i= i 



where 



for any Borel subset u C S n . By weak continuity of the integral Gauss curvature, 

/ f(af(N))da(N) = [ f(x)dp(x), (39) 
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which implies (37). 



On the other hand, for any pair (h, p) G A we have 

h(N) - p(a-/(N)) > }C(af(N), N) = ti(N) - P(af(N)) ViV G S" (40) 

Integrating the left hand side of this inequality against the measure da(N) and using the 
change of variable formula (39) we obtain 

/ h(N)da(N)- I p{af{N))da{N) = I h(N)da(N) - I p{x)dp{x) = Q[h, p]. 
Js n Js n Js n Js n 

Integrating the right hand side of (40) gives 

/ h(N)da(N)- I p(af(N))da(N) = I h{N)da{N) - I p{x)dp{x) = Q[h, p\. 
JS" Js n Js n Js n 

Thus Q[h,p] < Q[h,p] V(h,p) G A. This proves (36). QED. 

Theorem 8. The minimizing pair in Theorem 7 is unique up to rescaling with any constant 
c> 0. 

Proof. Let (hi,pi), i = 1,2, be two minimizing pairs which are not constant multiples of 
each other. Rescale, if needed, one of the pairs so that 

p 2 (x) < pi(ar) G S" 

with the equality holding for at least one x. It follows from (12) that h,2{N) < h\(N) ViV G 
S n . Assuming that p 2 ^ pi, it is clear that there exists N such that /^(^O < hi(N). Since 
the support function of a convex body is continuous, this inequality holds for some set of 
positive measure on S n . Then 

Q[h2,Pi] < Q[hi,pi]. 

This contradicts the minimality of (hi,pi). Therefore, h 2 (N) = hi(N) almost everywhere 
on S™. 

Let us show that h 2 (N) = hi(N) \/N G S™. Extend both functions from S n to the 
entire R n+1 as positively homogeneous functions of order 1 and denote these extensions by Hi 
and H 2 . The functions Hi and H 2 are locally Lipschitz and differentiate almost everywhere 
[34], section 1.5. Therefore, dH 2 — dH 1 = almost everywhere, that is, H 2 — Hi = const 
almost everywhere. Then H 2 — H 1 =0 almost everywhere. Since H 2 — H x is locally Lipschitz, 
we conclude that H 2 — Hi = everywhere. This implies that h 2 — hi = on S n and 
consequently pi = p 2 . QED. 

Remark 9. The convergence of minimizers P K to a minimizer of Q[h, p] combined with the 
uniqueness of the minimizer provides a constructive way for finding the solution to Aleksan- 
drov's problem numerically. 
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5 Aleksandrov's problem and optimal transport on S 



5.1 Connection with the problem of Monge 

Consider the problem of finding a map 9 E Q such that 

sup/" c(0- 1 (N),N)da(N)= [ c^^N), N)da(N), (41) 
6>e© Js n Js n 

where G is defined by (7)-(9) and c by (6) in the introduction. Note that 0^0. For example, 
the generalized Gauss map ap of the minimizer in Theorem 7 is in 0. The required properties 
of ap follow from known results in convexity theory (see, for example, [34], section 2.2) and 
the change of variable formula (39). 

The maximization problem (41) can be viewed as a variant of the celebrated Monge 
problem [26] . The original problem of Monge is formulated for two measurable sets U and V 
in Euclidean space R™ and two Borel measures m\j and my with equal total masses defined 
on U and V, respectively. It consists in finding a map To among all measurable maps {T} 
of U onto V such that m u (T~ 1 (B)) = m v (B) for all Borel sets B C V (that is, T pushes 
m\j forward to my) , and 

inf / l(x,T(x))dmu = / l(x,T (x))dmu, 
i T } Ju Ju 

where / is the Euclidean distance; see, for example, [17], Ch. VIII, section 4, or [36]. In 
our case the problem (41) is considered on S n and the cost function c(x, N) is nonlinear and 
allowed to assume infinite values. 

Theorem 10. The problem (4-1) admits a solution 6q and any such solution satisfies O^ 1 = 
a^ 1 almost everywhere on S" with respect to a (a — a.e.). Here, F G T n is the unique (up to 

rescaling) convex hypersurface defined by the unique (up to rescaling) minimizing pair (p, h) 
of the functional Q in Theorem 7. Furthermore, 

Q[h,p}= I c{af{N),N)da{N) = sup / c(0 _1 (AO, N)dcr(N). (42) 
Js™ fee is™ 

In addition, for any Borel set V C S™ for which fi{V) > we have 9 (V) = apiV). 

Proof. Let, as before, (h, p) be the pair in A defining F. For an arbitrary pair (h, p) G A 
and any 9 G we have for almost all iV G S n such that (0~ l (N), N) > 

h(N) - pie-^N)) > /C(e~\N), N). (43) 



16 



Integrating against dcr(N) and using (9), we obtain 

Q[h,p(9- l )]> [ JC(0- 1 (N),N)da(N) > I c(9~ l (N), N)da(N). (44) 

J{Nes« | (e- 1 (N),N)>o} is™ 

On the other hand, if h = h, p = p and 9 = ap then 

Ji(N) - p{af{N)) = t(af(N), N) ViV e S". (45) 
Since ap G O, the right hand side in (44) attains its supremum on 9o such that 

(aj\N),N) = (9 1 (N),N) cr-a.e. on S n . (46) 

Let us show that Oq 1 = a~p) a— a.e. on S n . It follows from (45) and (46) that 

m\N)) = P{uf{N)) <r-a.e. on S". (47) 

Let M be the union of the sets where either a^ 1 or 9$ l is multivalued and let N E S n \ M. 
By Theorem 8 the map ap is defined uniquely and then by (46) and (47) 

(p(9 o \N))0o\N),N) = (p( a ?(N)) a ?(N),N). 

This means that the points p(9 1 (N))9q 1 (N) and p(a^ 1 (A^))o;^ 1 (A^) lie on the same hyper- 
plane with normal iV supporting to F. Since F is convex, it contains the linear segment 
joining these two points and this linear segment is also contained in the same hyperplane. 
If ap 1 (N) 7^ 9^ l (N) then this segment does not reduce to a point and then ap 1 (N) is not 
single valued. This contradicts the choice of N. Therefore, 9^ = a^ 1 a— a.e. on S n . This 
and (45) imply that we have equalities in (44) when 9' 1 = a^ 1 . Consequently, (3) holds. 

Now we prove the last statement of the theorem. Let V C S™ be a Borel set and 
p{V) > 0. By (37) cr(ap(V)) = p{V) > and then af(ap(V)) = 9 \ap(V)) a- a.e. in 
ap(V). Therefore, 

V = Vnaf{ap{V)) = Vr\9 Q \ap{V)). 
Applying #o to both sides of the last equality, we obtain 9 (V) = ap(V). QED. 

Remark 11. A curious feature of the problem (41) is that it is, in fact, an intrinsic problem 
on the sphere S n in the sense that the class of admissible maps and the cost functional are 
expressed in terms of the logarithm of the cosine of the geodesic distance on S n . Though the 
condition (2) in Theorem 1 is formulated here for convenience as an extrinsic condition it is 
in fact intrinsic as the notion of polarity (= duality) in S n and (2) can be stated intrinsically 
[35]. Remarkably, the optimal map is (almost everywhere) the generalized Gauss map of a 
closed convex hypersurface in R n+1 . Thus, the solution of the optimization problem (41) on 
S n is a solution of the problem of embedding S™ into R™ +1 with prescribed integral Gauss 
curvature. 



17 



5.2 A representation of the generalized Gauss map 



For the optimal pair (p, h) and the corresponding convex hypersurface F G T n the equation 
(14) defines implicitly the generalized Gauss map a p. In this section we give an explict 
representation for this map in a form resembling the representation of optimal maps between 
subsets in Euclidean space with convex or concave cost functions [6], [7], [10], [11], [8]. 
However, in our case we give a representation for ap valid at all points of S™ and the 
spt/x is allowed to be atomic. This stands in contrast with usual assumptions under which 
representations of optimal maps are derived [11]. In order to state our result we need some 
definitions. 

Definition 12. Let f G C(S n ), x G S" and TS" o the tangent space to S™ at x . The set 

df(x ) = {ve TS n XQ | f(x){-v + f(x )x ,x) < f 2 (x ) \/x G S"} 
is called the subdifferential of f at x . 

The following proposition clarifies the geometric meaning of a subdifferential in the 
case of the radial function p of a convex hypersurface F G T n . 

Proposition 13. Let F G T n and p its radial function. For an arbitrary x G S™ denote 
by M(x) an outward normal to a hyperplane supporting to F at r(x) = p{x)x rescaled (if 
needed) so that (M(x),x) = p{x). Then v = —M(x) + p(x)x G dp(x). Conversely, for each 
v G dp(x) there exists a unique hyperplane supporting to F at r(x) with outward normal 
M(x) such that v = —M(x) + p(x)x. 

Proof. Fix some xq G S n and let P be a hyperplane supporting to F at ro = p(xo)xq. The 
set of supporting hyperplanes to F is not empty at any point of F and such P exists. Let M 
be the outward normal to P. Denote by Mq the projection of M onto the one-dimensional 
subspace generated by x and let Mg 1 " = M — Mq. Note that (M ,x ) > because the 
origin O is strictly inside the convex body bounded by F and M is an outward normal. Put 

M< = pMM l 
° (M Q ,x ) 

Then (Mq,x ) = p(x ). Since P is supporting to F at r we have 

(p(x)x, Mq) = p(x)(x, Mq T + Mq- 1 ) < (p(x )x , Mq) = p 2 (x ) Vx G S n . 

Thus, v := — Mq T G dp(x ). This also shows that for any x G S™ the projection on TS" of 
an appropriately rescaled outward normal to any supporting hyperplane to F at p(x)x gives 
a uniquely defined element in dp(x). 

Conversely, let v G dpo- Consider a hyperplane with normal vector M = — v +p(xo)x 
containing the point r . It follows immediately from the definition of v that this plane is 
supporting to F at r . QED. 
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Remark 14. At points of differentiability of F where there exists only one supporting hyper- 
plane the function p is differentiate and dp(x) = {gradp(rr)} with gradp(rr) G TS™ 7 where 
the gradient is computed with respect to the standard metric ofS n . By Rademacher's theo- 
rem almost all points on F (in the sense of Lebesgue measure on F) and consequently on S n 
(since it is a radial projection of F) are points of differentiability [34], notes to section 1.5. 

Corollary 15. Let F G T n and p its radial function. Then for each x G S n and each 
N(x) G cxp(x) there exists a unique v(x) G dp(x) such that 

N(x) = (48 ) 

y/V 2 (x) + p 2 (x) 



and 



<mA.v)= { -vW + pW% v (x)edp(x)S. (19) 
^Jv 2 (x) + p 2 {x) 



5.3 Connection with the problem of Kantorovich 



In this section we show that in the language of the optimal mass transport theory the 
minimization problem (18) can be viewed as the dual of the following primal problem. Denote 
by T(p, a) the set of joint Borel measures on S™ x S™ with marginals p and a, where p is as 
in Theorem 1 and a, as before, the standard Lebesgue measure on S n . Thus, any 7 G r(/i, 0) 
satisfies 

7[17, S n ] = p[U] and 7 [S n , U] = a[U] for any Borel U C S". (50) 
The primal problem in this setting is to determine a 7 G r(/i, a) such that 

C^} = sup C[ 7 ] (51) 

r( M ,<r) 

where C[ 7 ] was defined by (10). The duality relation between the problems (51) and (18) is 
given by the following 

Theorem 16. Let F(p,a) and C be as above. Then there exists a 7 satisfying (51) and 

Q[h,p[=C[% (52) 
where (h, p) is the optimal solution of problem (10). 

Proof. It is clear that the Gauss map of the optimal solution F = (h,p) in Theorem 10 
gives a measure 7 D G r(/i, a). In fact, any map from the set defined in section 5.1 gives a 
measure in T(p, a). To see this, note that if 6 G 6 then one can take je[U, V] = a[9(U) fl V] 
for any Borel sets U, V on S™. Then by (9) 

Che] = [ c(9- 1 (N),N)da(N). 
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Consequently, taking into account (42), we get 

sup C[i\ > sup I c(0- 1 (N),N)da(N) = [ c(aJ 1 (N),N)da(N) = Q[h,p\. (53) 
r(n,a) dee Js« Js™ 

We prove now the reverse inequality. Since c(x, N) < V(rr, N) G S n x S n , it is clear 
that in the maximization problem (51) it suffices to consider only 7 G T(/i, a) with support 

spt'y C {(x, N) G S" x S n I (x, N) > 0}. (54) 

Denote by r + (/i, a) the subset of measures in r(/i, a) satisfying (54). For any pair (h,p) 
from the set A of admissible functions (16) and any 7 G T + (fi, a) we obtain, taking into 
account (17), 

C[i]< [ [ [h(N)-p(xMdx,dN)= f h(Nh(S n ,dN)- [ P (xMdx,S n ) = Q[h,p\. 

(55) 

It follows from (53) and (55) that the supremum in (51) is attained on 7 corresponding to 
(h, p), that is, 7o = 7 and (52) holds. QED. 
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